Discovery of Three Distant, Cold Brown Dwarfs in the WFC3 
Infrared Spectroscopic Parallels Survey 



D. Masters 1 ' 2 , P. McCarthy 2 , A. J. Burgasser 3 , N. P. Hathi 2 , M. Malkan 4 , N. R. Ross 4 , B. 
Siana 1 , C. Scarlata 5 , A. Henry 6 , J. Colbert 7 , H. Atek 7 , M. Rafelski 8 , H. Teplitz 8 , A. 

Bunker 9 , A. Dressier 2 

ABSTRACT 

We present the discovery of three late type (>T4.5) brown dwarfs, including a 
probable Y dwarf, in the WFC3 Infrared Spectroscopic Parallels (WISP) Survey. 
We use the G141 grism spectra to determine the spectral types of the dwarfs 
and derive distance estimates based on a comparison with nearby T dwarfs with 
known parallaxes. These are the most distant spectroscopically confirmed T/Y 
dwarfs, with the farthest at an estimated distance of ~400 pc. We compare 
the number of cold dwarfs found in the WISP survey with simulations of the 
brown dwarf mass function. The number found is generally consistent with an 
initial stellar mass function dN/dM oc M~ a with a = 0.0-0.5, although the 
identification of a Y dwarf is somewhat surprising and may be indicative of either 
a flatter absolute magnitude/spectral type relation than previously reported or 
an upturn in the number of very late type brown dwarfs in the observed volume. 
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Introduction 



Understanding the nature and demographics of low-temperature brown dwarfs, which 
bridge the gap between dwarf stars and massive planets, is an important go al of modern as- 



trophysics. Since the discovery of the first methane dwarf G1229B in 1995 flNakajima et al. 



1995 



Y ( jCushing et al. 



Qppenheimer et al.l 119951 ). the number of known ultracool dwarfs of type T and now 



20111) has continued to grow, primarily due t o deep n ear-infrared surveys 



such as 2MASS fcirkpatrick et~al~1ll999l: IBurgasser et al1l2002l). SDSS (Strauss et a.. 



Chiu et al.ll2006l). UKIDSS (ILodieu et al. 



2007 



Chiu et al. 



2008 



Burningham et al. 



1999 



20101) 



CFBDS flDelorme et al.ll2008f ) and WISE flWright et all 



2010; 



Kirkpatrick et all 120111) Ad 



ditional dwarfs have been identified i n proper motio n surveys (e.g. iKirkpatrick et al.l 12010 



Lucas et al. 


2010; 


Scholz et al. 


2011; 


Liu et al. 


2011) 



finding brown dwarfs in the solar neighborhood; however, due to their faintness (a T8 dwarf 
has Mab ~ 17 in H band), probing the larger-scale galactic distribution of cold dwarfs has 
not been possible. 

Deep surveys in the near-infrared using the Wide Field Camera 3 (WFC3) on the HST 
enable the discovery of cold dwarfs at much larger distances. WFC3 parallel im aging surveys 
have identified cold dwarfs based on broad-band colors (e.g. iRyan et al.l 120 111 ), but the lack 
of spectroscopy means that the spectral types are uncertain. WFC3 near-infrared slitless 
grism spectroscopy, in contrast, permits unambiguous identification and classification of T 
dwarfs, whose sp ectra show pronounced atmospheric CH 4 and H 2 absorption bands in 
the near-infrared (IGeballe et al.l 12002c IBurgasser et al.ll2002l ). Moreover, the WFC3 spectra 
reach depths that are unachievable with ground-based surveys. 

Here we present three of the most distant cold brown dwarfs known, discovered in the 
WFC3 Infrared Spectroscopic Parallels (WISP) Survey. These discoveries illustrate that, 
by enabling spectroscopic identification of ultracool dwarfs at large distances, HST-WFC3 
grism spectroscopy can begin to probe the galactic spatial distributio n of the ultracool dwar f 
population and search for substellar members of the halo population ( IBurgasser et al.ll2003l ). 



2. WISP Survey Overview 



The WISP Survey ( lAtek et al.l |2010| ) has obtained slitless spectra over 0.8-1.7/im for 
more than 700 arcmin 2 of sky using the two infrared grisms installed on the IR channel of 
Wide Field Camera 3 (WFC3). While the primary goal of WISP is to measure the star 
formation rate over 0.5 < z < 2.5, the wide area and variety of galactic latitudes covered by 
WISP allow for the discovery of late-type dwarfs out to > 1 kpc. 
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WISP is a "pure parallel" program, meaning that the pointings are determined by 
other observing programs using either the Cosmic Origins Spectrograph (COS) or the Space 
Telescope Imaging Spectrograph (STIS). The WFC3 observations are taken in a parallel 
field determined by the offset of the instruments and the roll angle of the telescope. WISP is 
divided into a deep, narrow survey for parallel targets with more than four orbits of visibility 
and a shallow, wide survey for parallel targets with one to three orbits of visibility. Parallel 
observing targets are selected at galactic latitudes |6| > 20° with a preference for longer 
visibility times. Typical integration times for a 4-5 orbit target are ~5000 sec in G102 and 
~2000 sec in G141. The shallower WISP fields are restricted to the G141 grism, with a 
typical integration time for a 2-orbit target of ~4000 sec. Because the wide-shallow survey 
actually achieves deeper G141 integrations and covers a larger area than our deep-narrow 
survey, faint dwarfs are more likely to be found in the wide-shallow survey. 



3. Ultracool Dwarfs Discovered in the WISP Survey 

Brown dwarfs of type >T2 are easily recognizable in H-band spectroscopy due to the 
presence of two prominent, broad pseudo-emission features at roughly 1.26 and 1.58 /im 
interspersed with deep absorption bands of water and methane in their atmospheres. The 
155 G141 grism exposures (each 123" x 136") currently obtained by the ongoing WISP survey 
were systematically scanned by at least two observers to search for these objects, resulting 
in the identification of the three sources shown in Figure [TJ Their properties are summarized 
in Table 1. 



4. Spectral Classification 



We determined the spectral types of the objects by co mparing the WFC3 s pectra to 
69 T dwarf spectra from the Spe X Prism Spectr a l Libr aries (IBurgasser et al.ll2006l ) and two 
additional Y dwarf spectra from ICushing et al.l (120111 ). which were also obtained with the 
HST/WFC3 G141 grism using the same reduction software as the WISP survey spectra 
presented here. The WISP spectra and brown dwarf template spectra were compared using 
a x 2 statistic. The classifications and associated uncertainties were then found by computing 
the weighted mean and deviation from the mean with a weighting factor of e^ mm ^ x2 ^ x2 ^ 2 . A 
cutoff weight of >0.1 was used so that the results are not skewed by a large number of poor 
fits to early-type T dwarfs (a sampling problem that arises because the template sample 
contains disproportionately more early- type T dwarfs). 
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WISP 1232-0033 and WISP 0307-7243 are found to be type T7±0.6 and T4.5±0.4, 
respectively. WISP 1305-2538 is found to be of type T9.5+, and is probably a Y dwarf. The 
spectrum of the Y0 dwarf WISE 1405+5534 is shown overlaid in Figure [U in addition to the 
T9 dwarf UGPS 0722-0540. While neither template is perfect, the Y0 spectrum is a better 
fit. 



Distance Determination 



With accurate spectral types derived from the WFC3 spectra, the distances to the dwarfs 
can be found by comparing their measured broad-band magnitudes with those of nearby T 
dwarfs with known parallaxes. For this comparison, we converted the WFC3 F140W/F160W 
photometry for the WISP dwarfs into MKO J or H band with color corrections computed 
using the best-fitting T dwarf template spectra, applying the same weightings as for the 
spectral type determination. Given the inferred MKO J/H m agnitudes, the distan ces were 
found using the absolute magnitude/spectral type relations of iDupuy fc Liul ( 120121 ). 



WISP 1232-0033 (T7) is found to be at 270±60 pc and WISP 0307-7243 (T4.5) at 
400±60 pc. These are the most distant T dwarfs currently known, and have vertical 
scale heights of approximately 240 and 260 pc, respectively. WISP 1305-2538 is more 
difficult, as only one Y dwarf has reported a preliminary (and high ly uncertain) parallax 
(WISE J154151. 66-225025.2, vr = O'.'35+O'.'ll: fcrkpatrick etaDboilh and the classification 
scheme for these objects remains in progress ( Cushing et al. 201lh . We used the known 



parallaxes of the T9 UGPS 0722-0540 (0 / /246±0'. / 0033; Lucas et al. 2010) and the T9.5 CF- 
BDSIR J145829+101343B (0"340±0"0026; Liu et al. 2011) to infer a distance of 40-60 pc 
for this source, albeit with considerable uncertainty because of the classification. 



Table 1. Cold Brown Dwarfs Discovered in the WISP Survey 



Name 


RA (J2000) DEC (J2000) WFC3 magnitude MKO magnitude Spectral Type Distance Galactic Latitude 




(h m s) (° ' ") 


(AB) 


(AB) 




(PC) 


(deg) 


WISP 1305-2538 


13 05 25.51 -25 38 28.8 


23.01±0.04 b 


22.33±0.04 d 


T9.5+ 


~40-60 a 


+37.1 


WISP 1232-0033 


12 32 42.42 -00 33 06.7 


23.34±0.05 b 


22.65±0.06 d 


T7±0.6 


270±60 


+63.1 


WISP 0307-7243 


03 07 41.12 -72 43 57.5 


22.67±0.02 c 


22.36±0.02 e 


T4.5±0.4 


400 ± 60 


-40.9 



a The interval arises from classification uncertainty and the current lack of quality parallaxes for Y dwarfs. 

b F140W 

C F160W 

d J band 

e H band 
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6. Space Density 



The dwarfs discovered in the WISP survey provide the first opportunity to examine 
the mass function and spatial distribution of these late-type dwarfs in a deep field (see 
Pirzkal et al.ll2009l for an equivalent study of late M and L dwarfs). To assess what constraints 
our small sample makes on these statistics, we computed the expected suface densities of T 
dwarfs as a function of spectral type in the WISP survey fields as follows. 



1. V olume dens i ties for the cold dwarf subtypes are computed using the methods outlined 
Burgasserl ( 12004 . 120071 ) . We assume a power-law stellar mass function dN/dM ~ M~ 



in 



with a varying from -1 to 1, nor malized to agree with measured field values in the 
0.09-0.10 M range (0.0055 pc 



-3. 



Reid et al. 


1999; 


Chabrier 


2001) 



an age distribution consistent with a constant star formation rate over 10 Gyr, and 
computed the l umino sity distribution for the dwarfs using the evolutionary models of 



Burrows et al.l (119971) . Finally, we use empirical bolom etric correction/spe ctral type 



( jLiu et al.ll2010l ) and absolute magnitude/spectral type (IDupuy fc Liull2012l ) relations 
to convert the luminosity function to a spectral type surface density distribution. 

2. The expected surface density of dwarfs per WISP 123" x 136" pointing is then found 
by computing the limiting magnitude, and thus the maximum distance, for which 
elcb subtype would be spectroscopically detected in a given field. The online WFC3 
exposure time calculator is used for this analysis. Corrections to the derived limiting 
magnitudes are applied to account for Malmquist bias. We assumed as a first order 
approximation that 10% of each pointing is lost due to overlap, edge effects, etc. as 
estimated from the galaxy extractions. 

3. The effective volumes for each field/subtype are found and multiplied by the space 
densities per subtype. The results for all the WISP fields are added, yielding a surface 
density for each subtype for the entire survey (Figure 2). 



The expected surface density of each spectral type across the entire WISP survey is 
summarized in Figure 2 for different assumed power-law slopes in the low-mass stellar mass 
function. The number of dwarfs found is roughly consistent with expectations; in particular, 
the fact that we find no early t ype (<T4) dwar fs, despite their higher intrinsic luminosities, 
is consistent with their rarity (IBurgasserl 120071 ) . The fact that we identified three T dwarfs 
is consistent with other field mass function estimates based on wide, shallow surveys of T 



1 We require S/N=10 in this analysis, an approximation based on simple tests with the WFC3 spectra. 
However, changing the S/N requirement can significantly influence the expected surface density of dwarfs. 
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dwarfs (e.g. iMetchev et al.l l2008r iBurningham et al.ll2010l ). However, finding a T9.5/Y0 
dwarf is somewhat surprising due to their faintness. Possible explanations for this include 
(1) a flatter absolute magnitude /spectral typ e relation than preliminary estimates indicate 
( ICushing et al.l 1201 ll ; iKirkpatrick et al.l 1201 if ). (2) a sudden upturn in the space density of 
very late-type dw arfs, possibly associated with a nearby free-floating planetary population 
( jSumi et al.ll201ll ) or a larger proportion of old halo brown d warfs at higher scal e heights (for 
this simulation, a 1:400 halo fraction was assumed based on lDigby et al.l 120031 ). (3) a fortu- 
itous discovery, or (4) some combination of the above. We note that WISP 1305-2538 is 4.2 



from 1 


;he very low-mass t 


2005 




Stumof et al. 


2008) 



Stumpf et al.ll2008f ). With a 5000 AU projected separation, this matching is unlikely 



( IDhital et al.ll2010l ). but second-epoch imaging would allow us to affirm or refute the exis- 
tence of a widely-separated, very low-mass quadruple. Additionally, follow-up proper motion 
measurements of all three new dwarfs would help establish whether they are galactic halo 
members. 



7. Summary 

We have presented three distant T/Y dwarfs spectroscopically confirmed by the WISP 
Survey. These discoveries show the power of HST-WFC3 slitless grism spectroscopy to si- 
multaneously find these objects at much larger distances than has been possible with ground- 
based surveys and yield accurate spectral types, making surveys such as WISP powerful new 
probes of the population of cold brown dwarfs in the galaxy. 



We thank an anonymous referee for helpful comments that improved this letter. This 
research has benefitted from the SpeX Prism Spectral Libraries, maintained by Adam Bur- 
gasser at http:/ / www.browndwarfs .org/ spexprismj We thank Michael Cushing for providing 
Y0 spectra from the WISE survey. 



This work has been supported by the Carnegie Observatories Graduate Research Fel- 
lowship. 
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Fig. 1. — HST-WFC3 G141 1D/2D spectra of the three cold brown dwarfs found in the 
WISP survey, with best-fit template spectra overlaid. Prominent pseudo-emission features 
at roughly 1.26 and 1.58 fim are characteristic of T dwarfs and make them relatively easy 
to identify in the 2D grism images. While WISP 1232-0033 and WISP 0307-7243 (middle, 
bottom) are in nearly perfect agreement with the best-fit templates, WISP 1305-2538 (top) 
is not perfectly fit by either a T9 or Y0 template, although Y0 fits better. The template 
spectra are from the SpeX prism library (IBurgasser et al.l (120101 ) and references therein). 
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Fig. 2. — Left: Expected cumulative surface density of late-type dwarfs across all WISP 
fields for different values of a in the power-law stellar mass function dN/dM ~ M~ a , 
from the simulations described in §6. The total number found (three, black line) is generally 
consistent with c*=0.0-0.5. Right: The expected surface density of each spectral type across 
the entire WISP survey. The spectral types of the three dwarfs found are indicated with 
downward arrows. The discovery of a T9.5/Y0 dwarf is somewhat surprising; see §6 for a 
discussion of possible explanations. 
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